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Abstract 

We decorate the square lattice with two species of polygons under the constraint 
that every lattice edge is covered by only one polygon and every vertex is visited by 
both types of polygons. We end up with a 24 vertex model which is known in the 
literature as the fully packed double loop model (FPL^). In the particular case in 
which the fugacities of the polygons are the same, the model admits an exact solu- 
tion. The solution is obtained using coordinate Bethe ansatz and provides a closed 
expression for the free energy. In particular we find the free energy of the four 
colorings model and the double Hamiltonian walk and recover the known entropy 
of the Ice model. When both fugacities are set equal to two the model undergoes 
an infinite order phase transition. 
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1 Introduction 



Loop models appear in a variety of contexts: as diagrammatic expansion of the Potts 
and the 0(n) models |2H1, and as models to describe polymers in several phases. The 
conformal properties were investigated using the Coulomb gas mapping ^T]. An exact 
expression for the free energy, in the thermodynamic limit, was found using the coordinate 
Bethe ansatz method and solving linear integral equations fHl ^] • Since then, new loop 
models have been studied on the square lattice [211 C] , on decorated lattices E2] , on 
a random lattice and in three dimensions The resolution techniques have been 
improved both on the numerical side, introducing the connectivity basis fI^ |23 and 
on the analytical front adding the seam fSl 02] , making use of the algebraic Bethe ansatz 
PBj . mapping to an interface model [30] and exploiting conformal field theory pU] . 

Here we are interested in the fully packed double loop model (FPL^) on the square 
lattice. It is defined by filling the square lattice with loops in two colors, drawn along the 
lattice edges in such a way that every bond is covered by one loop and every site is visited 
by both types of loops. It was introduced for the first time in as an alternative repre- 
sentation of the four colorings model on the square lattice. It is a natural generalization 
of the fully packed one loop model on the honeycomb lattice j23 1201 1201 Ej , which is the 
loop representation of the three colorings model on the hexagonal lattice Jlj . 

Jacobsen and Kondev mapped the FPL^ model onto a height model and postulate that 
the long wavelength behavior is correctly described by a Liouville field theory ^Uj. This 
provides a geometrical view of conformal invariance in two dimensional critical phenomena 
and a method to access the critical properties of loop models exactly. In particular, they 
succeed in calculating the central charge and the critical exponents ^IHI- 

The FPL^ model exhibits a rich phase diagram and provides a representation for 
previously studied models: Ice model [lOj, four colorings model [01, dimer loop model 
|3Uj . double Hamiltonian walk [lOJ, compact polymers |221- A generalization of the FPL^ 
model, which offers a unifying picture of the compact, dense and dilute phases of polymers, 
was proposed in |4j by relaxing the full packing constraint. Another possible generalization 
corresponds to the Flory model of polymer melting p|. Recently the FPL^ model has been 
coupled to two dimensional quantum gravity in order to study meanders |31j . 

Here we find an exact expression for the free energy, in the thermodynamic limit, 
when the loop fugacities are the same. The solution relies on coordinate Bethe ansatz. 
Unfortunately the compact polymer does not belong to the solvable line. The configura- 
tional entropy for the double Hamiltonian walk |10[ and the four colorings model [0] is 
calculated. We also confirm a conjecture on the average length of loops [2|. 

The present work is organized as follows. In section 2 we map the FPL^ model into 
a 24 vertex model on the square lattice. The transfer matrix is set up in section 3 and 
diagonalized by means of coordinate Bethe ansatz (BA) in section 4. The free energy, in 
the thermodynamic limit, is calculated in section 5. 
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2 The model 



We study a 24 vertex model on the square lattice. Every edge is decorated with one of two 
types of arrows, black or grey ones. The admitted configurations of arrows do not contain 
sources and sinks, so that for both colors there is one arrow going in and one going out of 
each vertex. The 24 possible vertices may be divided into 6 basic types depicted in Fig. HI 
the others being related by 7r/2 rotations. We can interpret each vertex as a couple of 
lines turning clockwise or anticlockwise. In this way the vertex model is mapped into an 
oriented loop model. We are interested in the partition function of the unoriented loop 
model defined as: 

where and Ng are the respective numbers of black and grey loops. We assign the 
fugacity rib for black loops and ng for grey loops. The sum is over all the allowed loop 
configurations. The fugacities of the loops are related to the weights of the vertex model 
by: 

= exp(i7refe/4) r/ = exp(i vr 6^/4) nb = 2 cosireb ng = 2cosTTeg (2) 

This is the two- flavor fully packed loop model (FPL^) on the square lattice investigated by 
Jacobsen and Kondev Ij. We will see that when the vertex weights are the same u = rj, 
the partition function of FPL^ (Q) can be computed exactly, in the thermodynamic limit, 
making a coordinate Bethe ansatz. For particular values of the fugacity we recover the 
partition function of previously studied models. In particular (n^, Ug) = (2, 2) corresponds 
to the four colorings model on the square lattice For {nh,ng) = (1, 1) we get the Ice 
model solved by Lieb jT3]- The single Hamiltonian walk is obtained for {nb,ng) = (0, 1) 
Finally {nb,ng) = (0,0) is the double Hamiltonian walk for which the entropy has 
been correctly conjectured in [T], as we will show. 




l/coT] {o/Tj COT) ri/ca 1 1 



Figure 1: The 6 basic types of vertices of the 24 vertex model and the corresponding 
weights. When uj = r] the model admits an exact solution. 
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3 The transfer matrix 



Here we set up a transfer matrix for the model following ^I^ I14j. Consider a square lattice 
on a cylinder, made up by LM vertices and connect the vertices by vertical and horizontal 
bonds. Starting from the bottom we have a circular row of L vertical bonds followed by a 
row of L horizontal bonds and so on alternately. Denote by |\E') the possible configuration 
of a row of L vertical bonds, that is to say a precise assignment of arrows on each of 
the vertical bonds. There are 4^ possible choices. Let \^) and |\E'') be the states of two 
consecutive rows of vertical bonds. The intervening row of horizontal edges is decorated so 
that the L resulting vertices are of the types shown in Fig.^ To each allowed configuration 
of vertices we assign a weight given by the product of the weight of each single vertex. 
Then we sum over all the possible configurations of horizontal bonds compatible with the 
vertical bonds and obtain the total weight denoted by (\E''|T|\I/). In this way we build 
up the transfer matrix T for the model. The physics of the model is encoded in the 
eigenvalues and the eigenvectors of the transfer matrix. If periodic boundary conditions 
are imposed in the vertical direction, the partition function is simply the trace of the 
transfer matrix raised to the power M: 

Z = TiT^ (3) 

In the thermodynamic limit the free energy f^o follows from the largest eigenvalue A(L) 
of the 4^-square matrix T by: 

/oo= lim ylogA(L) (4) 

If we interpret the vertical axis of the cylinder as time (running downward, for conve- 
nience) and the horizontal extend as space we end up with the time evolution of a one 
dimensional system. 

In the following we are going to study the possible evolution of a generic state. Impose 
periodic boundary conditions on the horizontal direction. It is convenient to start from 
the state which consists of only gray arrows pointing up (pseudo- vacuum) . This can 
match only with a successive row of gray arrows up. There are two possibilities for the 
intervening row of horizontal bonds: they must be black arrows, either running all to 
the right or all to the left. The total weight is two. The next step is to replace some of 
the gray arrows up by black arrows pointing up. We call a black arrow up an ordinary 
particle. In this case the upper row can evolve in two different lower rows as illustrated by 
Fig. 121 all the particles move either one step to the right (Fig. |21 a) or one step to the left 
(Fig. 121 b). The product of the weight of each single vertex equals one, because vertices 
with reciprocal weight come in pairs, so that the total weight of the transition is set to 
one. More complex patterns arise allowing for black arrows pointing down. We refer to 
a black arrow down as a special particle. Now, besides the global shift, there are other 
possible evolutions of which an example is depicted in Fig.|21c. The special particle can be 
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connected via a half-loop to the left or right neighboring ordinary particle. Then in order 
to fulfill the boundary constraint a black half-loop has to be inserted in the lower row. 
This half-loop can turn either clockwise or anticlockwise and must be located between 
any two particles other than the two already connected via the upper row. Thus we end 
up with one half-loop in the upper row and one half-loop in the lower row. Particles shift 
one step to the right or to the left according to their position relative to the half-loops. 
The total weight is completely determined by the orientation of the two half-loops. A 
clockwise half-loop contributes ixP' and an anticlockwise one tu"^. To complete the list of 
allowed configurations we have to examine the states containing grey arrows down. We 
will interpret the grey arrow down as a bound state between two black arrows (see Fig.Elb 
and Fig. E}. We leave the details to the next section. 

A fundamental ingredient, in seeking for an exact solution, is the existence of con- 
served quantities. Besides the conservation of black and grey flux, there is an additional 
conservation law, which is less evident. We saw that at every application of the transfer 
matrix, both ordinary and special particles, shift one step. Therefore the lattice can be 
divided into two sublattices such that the number of particles on each sublattice remains 
constant from row to row. As a result the transfer matrix splits in a series of diagonal 
blocks relating states in the same sector. Each sector is completely fixed by specifying 
the number of odd, even and special particles. We are going to diagonalize this matrix 
making an appropriate ansatz on the structure of the eigenvectors. 

Before concluding this section we make some remarks concerning the relation between 
the loop and the 24 vertex model. We assign the weights exp(±i vr e^^g) to each oriented 
loop. The orientation fixes the sign and the color selects the phase or . The partition 
sum of the oriented FPL^ model may be cast in the form: 

Z = y^[exp(i7refc) + exp(-i7re;,)]^''[exp(i7reg) + exp{-in Cg)]^" (5) 

where the sum runs over all the allowed unoriented loop configurations. Relating the 
unoriented loop fugacities to the oriented loop phases by = 2 cos neb^g, we recover 
the partition function of the unoriented loop model (H)). The advantage of splitting the 
loop fugacity among the two possible orientations, is that in this way the loop fugacity 
can be distributed along all the vertices visited by the loop. This is achieved assigning 
the weight exp(±i tt e{,,g/4) to every vertex visited by the oriented loop, with sign + for 
a clockwise turn and — for the anticlockwise one. Since for every oriented loop, on the 
square lattice, the difference between clockwise and anticlockwise turns is ±4, multiplying 
over all the local vertex weights we recover the correct global loop weight. This should 
clarify relations ^ and the choices made for the vertex weights in Fig. ^ It is important 
to notice that the previous argument does not work for a loop winding round the cylinder 
because in that case the difference between clockwise and anticlockwise turns is zero. In 
order to give the correct weight to non-contractable loops a seam [THl UHl 1201 is placed 
between two generic vertices. This is well depicted in Fig. |21 The total weight of a generic 
transition will factorize in the weight of each single vertex and an extra factor which 
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accounts of the seam. The additional factor is set equal to a (a^^) when the seam cuts a 
right (left) pointing arrow. For convenience we introduce a phase a setting a = exp(i7ra). 
In terms of the phases the partition function of the unoriented FPL^ model on the cylinder 
reads: 

Z = ^(2 cos 7ref,)^H2 cos TT 63)^^(2 cos Tra)^'' (6) 

where Na is the total (black+gray) number of uncontractable loops. Notice that, in order 
to simphfy the calculation, we do not distinguish between gray and black uncontractable 
loops. 

4 The Bethe ansatz 

Since the grey arrow down serves as bound state between two black arrows, we restrict, 
for the moment, to the sector containing only ordinary and special particles. Denote the 
generic state of a row of L vertical bonds by |x, r). The vector x = (xi, . . . ,xn) gives 
the position of all the particles while the vector r = (ri, . . . ,rm) specifies that those 
at position Xr-^ , • • • , are special ones. The eigenvalue equation for the transfer matrix 
reads: 

T J2 ^(x> r) |x, r) = A ^ ^(x, r) |x, r) (7) 

|x,r) |x,r> 

where the eigenstate of eigenvalue A is a linear combination, with coefficients \l'(x, r), of 
the base |x, r). In the coefficients the eigenvalue equation assumes the form: 

J2 (x, r|T|x', r')vl>(x', r') = Avl/(x, r) (8) 

|x',r') 

Here we examine the sector with A^ particles in total, (A^ — 1 ordinary + 1 special), located 
at position xi, . . . ,xn, where Xj+i — Xj > 2, and suppose that the jth particle is a special 
one. Moreover, we assume that the seam is placed between two ordinary particles. In 
order to solve the eigenvalue problem we make the following Ansatz on the form of the 
wave function: 

N 

<i/j{xi, ...,xn)=j2 ^pi-Piv^j(p) n (9) 

p 1=1 

where the sum runs over all permutations p = {pi, . . . ,Pn) of the numbers 1, . . . , A^. The 
goal is to chose Ap, ipj{p) and the complex numbers Zp^, . . . , Zp^^ so that (jSI) is satisfied. At 
this stage the structure of the coefficient ipj{p), the behavior of Ap^ . p^^ under the exchange 
of two consecutive indices, and the dependence of the eigenvalue A on , . . . , Zp^ , is still 
unknown. In the following we are going to determine it. 

Inserting the ansatz Q into the eigenvalue equation (jHJ, we see that a sufficient 
condition for the validity of ^ is that the coefficient ipj{p) obeys, for any permutation 
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Figure 2: Four ordinary particles (black arrows up) and one special particle (black arrow 
down). Dashed line indicates the position of the seam. The special particle, in the upper 
row, is the second one from the left {j = 2). The possible evolutions of the upper row 
are: a) Shift to the right, b) Shift to the left, c) Appearance of the black half-loops. In 
the lower row the special particle is the fourth one (m + 1 = 4). 

p, the following evolution equation: 

A^j = a zi.,N i'j-i H ^ — ■ipj+i 

j+N~2 ^ 

+ O;^ (t^^^am^m +t^^/3m+l^m+l) — (10) 

m=j ■' 
i+7V-l 

+ 1^"^ (u;~^am^m +^^^/3m+l^m+l)^^^^^^^ 

7Tt=J + l ■' 

Since we will be concerned mainly with one permutation p, we introduce the more compact 
notation Zi = Zp^ and Zi.k = Y['j=i ^pj^ together with the cyclic condition Zj = zj+n- Let us 
examine in details all the terms. The first term in the rhs of (fTUI) describes the process in 
which all the particles shift one step to the right (Fig.|21a). Notice that the special particle 
is the jth one in the upper row and the (j — l)th in the lower row. The weight of the process 
is a because the arrow crosses the seam in the right direction. The product Zi^jy accounts 
for the global shift to the right. The second term describes the left translation (Fig. |21b). 
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Let us clarify the two sums. Start with the first one (Fig. |21c). The special particle in the 
upper row is connected via a clockwise half-loop to the left neighboring ordinary particle. 
This is taken into account by the coefficient uP' in front of the sum. The half-loop in the 
lower row must be inserted between any two particles except the ones already connected 
by a half-loop in the upper row. Therefore the sum runs from m = j to m = j + N — 2. 
Moreover the particles distributed between the jth and the mth particle undergo a shift 
to the right while those positioned between m + 1 and j ' + — 1 move one step to the 
left. These translations are taken into account by the products zj^m/ Zm+i-.N+j-i- Since 
the half-loop in the lower row may turn clockwise or anticlockwise, each term in the sum 
splits in two parts: u'"^ ipm for the anticlockwise half-loop and uo"^ ipm+i for the clockwise 
one. The same argument applies to the second sum. The coefficients and Pm+i ensure 
that loops which wrap around the cylinder pick up the correct weight. They are set equal 
to a {a~^) when an arrow, of any of the two colors, in the intervening row of horizontal 
edges, crosses the seam in the right (left) direction. Their values depend on the relative 
position of the special particle in the upper row, the seam, and the half-loop in the lower 
row. Following Fig.|21c we see that when the lower half- loop is placed between the special 
particle j and the seam, the intervening arrow, cutted by the seam, points to the left, so 
that am = Pm+i = When it is positioned on the right of the seam = Pm+i = ci. 
In the special case in which the seam cuts the lower half-loop am = and f3m+i = d- 

Our aim is the computation of the coefficient ipj and the eigenvalue A. Rearranging 
the terms, equation ()10j) may be expressed in a more compact form as: 

X 

A -ipj = -a zi,N V'j 1 + [zj^ ^^'^ + Qj (11) 



where the quantity Qj is defined as: 



Qj= (u;"^am^m + V'm+i) ^-^^ — (12) 

m=j Zm+l:N+j-l 

The term Qj generates all the oriented black half-loops in the lower row, that can be 
accommodated between two consecutive particles. An alternative definition of Qj comes 
from equation (|TT|l : 

A + azi.NUj'^ + a-'^ z'l^u-"^ 
Qj = -1 2 ^^3) 

Zj ^ + Zj UJ^ 

From the original definition of Qj p2|) it follows that: 

- Qj = (^^~^^j + ^^^j+i) (azi;N — ] 

Substituting in the Ihs of (fT^ the alternative expression for Qj given by (fT^: 

(A + a zi,N uj^ + a"^ 2"^ uj'^) ( — r-^^^i ^j+i o] 2 2 



(14) 
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Figure 3: a) The two particles are at a distance of two, no black loop can be inserted in 
between, b) Particles on adjacent sites force a gray half- loop to be inserted. 



(a; ^pj +u i/jj+i) a Zi.,n 



(15) 



Eq. (|T5|) can be read as a recursive relation for ipj. It simplifies defining the eigenvalue 
via a new variable /i: 



A 



a Zi-i\j- 



{u' - /x) + 



azi-N 



(16) 



Replacing the eigenvalue A in Eq. (fT3|) by the definition ()16|) and dropping the factor 
which contains azi-^-, the recursive relation for ipj reads: 



io ) z : 



+ z: 



■i+i 



'i+i 



(17) 



Notice that relation ()17j) does not involve the parameter a. To see how the parameter a 
enters in the wave function we have to inspect the case in which there are no ordinary 
particles between the seam and the special particle. Notice that in this special case the 
seam may cut the upper black half-loop. We skip the calculation and give the final result 
for the structure of the coefficient 



-2 I , ,2 ^2 
+ UJ Z^ 



m—1 



n 



jj + uj + zl 



(/. - u;^) zl - 1 l\ (/i - uj^) zl 



18) 



where s = 1 if the special particle is located left of the seam, and s = if it is positioned 
on the right side. 

The formalism developed up to now, works for sparsely distributed particles. When 
the condition Xj+i — Xi > 2 does not hold, several restriction on the amplitude Ap are 
necessary. Suppose that two generic particles are located at position x — 1 and a; + 1 (see 
Fig. El a), then a half- loop of black arrows can not be inserted between them in the lower 
row. Thus, this event must be omitted from the generic list of possibilities, in which the 
particles have both came from afar. We demand the weight of the corresponding term. 
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automatically generated by the Ihs of the eigenvalue equation when the particles are 
sparse, to be zero in this special case: 

A...jk... (^"^ ttm V^m + UJ^ Pm+l V^m+l) = (19) 

where the sum is over the interchange j ^ k of the indices and the shorthand notation 
A...jk... = ^...pjPfc... has been used. The way in which the momenta Zj and Zk enter in 
the coefficients ipm and ipm+i can be read from the structure of (fTHj). Inserting (fTHj) in 
Eq. (fT^ and working out we see that only terms proportional to fi survive. The resulting 
requirement for the scattering factor is: 

^ 1 + (^^ + uj-')z] + z]zl 
A.k,... l + {co^ + co-^)zl + zpi ^ ^ 

Let us make the following important observation: the transfer matrix acts on a generic 
state shifting each particle one step to the right or to the left. Thus we can group particles 
in two sets according to their position: Odd and Even particles. The two families do not 
mix. This means that the scattering relation ()20j) holds for particles belonging to the 
same family. 

In the following we elaborate our formalism in order to incorporate states containing 
grey arrows pointing down. Start by noting that black arrows positioned on adjacent 
sites generate a half-loop of grey arrows (see Fig.Elb). This observation leads to treat the 
grey arrow down as bound state of two generic particles. Consider the states with N — 2 
ordinary particles (black arrows up) and 1 bound state (grey arrow down). We make the 
following natural Ansatz on the form of the wave function: 

$^(xi, . . . , . . . XN-i) = J2 \....pnMp) • • • (^p™ • • • z;^-'] (21) 

p 

where the index m indicates that the bound state is preceded by m — 1 ordinary particles. 
The bound state is located at position Xm- Notice that the wave functions (j?T|) and © 
have essentially the same structure apart from the coefficient 0m which is still unknown 
at this stage. In order to compute it we work out the eigenvalue equation (jSJ with the 
help of the ansatz ()21|) and for the particular case in which the initial state (rhs of 
Eq. ISI) contains two adjacent particles. We see that if we require that: 

A...jk..XuJ~^i^m + Uj'^i'm+l)Zj = ^ A,„jk...4>m{v'^ + V^ZjZk) (22) 

then the generic list of possible evolutions (see rhs of Eq. inUl automatically incorporates 
the gray half-loop. The coefficient t] accounts for the two possible orientations of the gray 
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Figure 4: Four ordinary particles (black arrows up) and one bound state (gray arrow 
down). The bound state is the second particle from the left (m = 2). The possible 
evolutions of the upper row are: a) The bound state go straight and the ordinary particles 
shift to the right, b) The bound state go straight and the ordinary particles shift to the 
left, c) The bound state forms a small gray half-loop and evolves splitting into two black 
arrows (the jth and kth particles). Periodic boundary conditions force a black-half loop 
to be inserted in the lower row. 



half-loop. The next step is to study the evolution of the bound state. The analog of 
Eq. (fTn|l for the bound state reads: 



ZjZk zi;N a, 
— \ 



(23) 



+ 



2\ m+N-1 

-2 I ^ \ Zm+2:n / -2 ; i 2 ^ ; \ 

Tj ^k + — ] 2^ Z anWn + UJ (3n+lVn+l) 



n=m+ 



^ Zn+l:m+N-l 



The bound state can move straight or form a half-loop as shown in Fig. Now we have 
to check the consistency of Eq. (|^. For that eliminate 0^ in Eq. (PHj) using relation 
(j22|) . In order to eliminate A, plug in the value of Kipj given by the rhs of the evolution 
equation (jl(J|) . At the end we get an expression relating ipj for different values of j. In 
the Appendix we will show that such expression is consistent if we choose the same value 
for the vertex weights a; = 77 (in term of phases = e^) and demand that the scattering 
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amplitude is symmetric for the interchange of particles belonging to different families: 



A...jk... = A,„kj... (24) 

Notice that the new scattering relation (|24p is not in contradiction with the previous 
relation ()2()|1 since the new one involves particles of different families. 

We went further and studied the sector with two special particles. We skip the calcula- 
tion and present the outcome. Unwanted terms arise in the expression for the eigenvalue. 
In order to cancel them we introduce another amplitude B^, in the wave function, which 
regulates the behavior of the special particles. We are forced to impose the following 
scattering relation: 

B...ij... _ jJ^ijJ^j + uj'^l^j - cuVi , . 

B...ji... fiifij +UJ ^fii- uj^fij 

Now we have all the ingredients to derive the Bethe ansatz equations. The final ansatz 
on the form of the wave function is: 

N 

^{x,,...,XNK...,rJ = J^^p^qJI-^p; (26) 

p,q j=i 



m 



n 



1=1 



in 



1 



The wave function has to satisfy periodic boundary conditions. Suppose that the first 
particle is an ordinary one (ri ^ 1): 

^(xi, . . .,XN\ri, ...,rm) = ^(X2, ...,Xn,Xi + L\ri - 1, . . . - 1) (27) 

From the Ansatz ()26p. it follows that a sufficient condition for the validity of (j^7j) . is that 
the complex variable Zp-^ fulfills, for every permutation p, the relation: 



^P...PN 11 ^^^^ _ ^4)^2^ _ 1 - ^P2-PNPi Si 



We make the additional assumption that the system has an even length L, this means that 
after having imposed periodic boundary conditions to a generic particle the sublattice to 
which the particle belongs does not change. We associate momenta zj {j = 1, . . . ,nz) 
with the even sublattice and i/j (j = 1, . . . ,ny) with the odd sublattice. By eliminating 
the amplitudes in using relations and ((211), we get the first two families of BA 
equations: 

+ + 4 - l + {u- + u;--)z] + zpi _ , 
11 (/.. - c.^).| - 1 n l + iuj^ + u-^)zl + z]zl 
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(30) 



Similarly, assuming that the first particle is a special one, we find an equation for the 
variable /i^^: 



Pk 



Making use of relation ()25p we get the third family of BA equations: 



(31) 



(32) 



The factor arises due to the fact that the special particle moves from the left to the 
right side of the seam when periodic boundary conditions are imposed. 
The eigenvalue expression that generalizes (fT?)|) has the form: 



A = (-1)"^ 



Thy Tin 



"'Z --y ■"H' -I ii'Z -I -"y -I ■■"H- -1 



j=l k=l m=l 



; = 1 k=l m=l 



(33) 



5 The free energy 

In order to find an exact solution of the model in the thermodynamic limit we need 
difference kernels ^H] • First simplify a bit notation (0) introducing the phase 6 = Cb = Cg 
and reminding that a = exp(i vr a). Then introduce a new set of variables uj, Vj,Wj related 
to Zj,yj,fij by: 



2 _ sin^(l + Mj i) 
~ sinf (1 -Mji) 



and define the function: 



2 _ sinf (1 + fji) 
~ sinf (l-t;,i) 



n f , sin^(c + xi 
Sc{x, e) = 



2 i sin TT 9 
1 — exp(7r6'wj) 



sm —{c — x\) 



The BAE in the new variables are: 

Uw riu 



i=l 



k=l 



(34) 
(35) 

(36) 
(37) 



k=l 
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exp(-2i7ra) ]^S'i(w/ - Mj, 6^) J]^ S'i(w« - Vk,9) J]^ -S2{wm - wi, 9) = -1 (38) 
j=i k=i 

The eigenvalue is the sum of two terms: 



in=l 



A 



riu riv 

exp(i vr a) JJ Si{uj, Qfl'^ JJ Sxivk-, 



77 sm— (lymi - 2) 



n 



k=l 



sm ^Wm 1 
m=l 2 



(39) 



+ expi-ina)llS,iu,,e)-'/^l[S^ivk, 

j=i k=i 



-1/2 -Q sinf (u;^i + 2) 



sm -TrWm 1 
m=l 2 



Following the argument of JH] we see that the dominant eigenvalue lies in the sector 
labeled by (n„ = Uy = Uw = L/2). For a = we get the largest eigenvalue when the 
roots are symmetrically distributed on the real axis in such a way that no holes appear 
between two consecutive roots (close packing of the roots). Switching on the twist a, the 
root distributions start to drift. We verify numerically that: 



j=i k=i 



a/2 



exp(i7raj 



(40) 



Which allows for an eigenvalue expression with just one term: 



m=l 



1/2 



(41) 



And the a dependence has been encapsulated in the factor A: 



an TT ^ coth 



m=l 



71 9Wr. 



(42) 



We can go further with the simplification taking the logarithm. For that purpose define 
the function: 

(l),ix,e) = i\ogS,ix,e) (43) 

And introduce the counting functions Zl{x,6) as follows: 



-01 (m, 9) + ^ 01 (Wm - M, d 

m=l k=l 



)i(w, 0) + "^ (pliwm -v,9)-'^ (piivk - V, 9) 



(44) 
(45) 



m=l 



k=l 
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TTiit ^^1) 



j=l k=l m=l 

Notice that the counting functions defined in this way have positive derivatives. The next 
step is to define the root density function via the derivative of the counting functions in 
the thermodynamic hmit: 

/ N 1 1- d Zl(x) 

2 TX L^oo dx L 

Numerically we find that approaching the thermodynamic limit the roots spread over the 
real axis, from — oo to +oo, thus the integral equations for the root density functions are: 

1 I 

Pu{u, 6) = -—(p'liu, 0) + — [4>[{x - u, 9)p^{x, 6) - (t)'^{x - u, 9)p^{x, 9)] dx (48) 
47r ivr 

1 1 /""^ 

pu,{w, 6) = — [(f)[{x - w, 9) {puix, 9) + p^{x, 9)) - 02(x - w, 9)p^{x, 9)] dx (50) 

The system can be solved using Fourier transform. The solutions for the root density 
functions are: 



Pw[X) 



^ sech (^^) p„(x) = py{x) = ch (^^) sech (^^) (51) 

It is remarkable to notice that they do not depend on the parameter 9. There is still 
9 dependence in the expression for the largest eigenvalue ()4H) . The free energy in the 
thermodynamic limit has the following integral representation: 

fooi9) = — / log cos^ n9 + sin^ n 9 coth^ sech [ — ] dx (52) 

16 y_oo L 2 J V 4 / 

Notice that the a dependence (j42j) drops out in the thermodynamic limit. Applying 
Parseval formula we can drop the logarithm in the integrand and write the free energy as 
a Laplace transform: 

. /""^ e-Psech(2p^) sinh(p^) ,1, 3, , , , 

For particular values of 9 the integral can be solved [T^ . 
Two mutually excluding Hamiltonian walks ^: 

/oo(l/2) = llog2 (54) 
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The Ice model [T^ : 

/oo(l/3) = ^log^ (55) 
Four colorings on the square lattice [01 Cj: 

foo{0 0) = -^log2 - ^ logvr + 21ogr(l/4) (56) 

Notice that in this last particular limit expression (j35|l becomes rational. These BA 
equations were derived for the first time in [E] and look very much like those derived in 
P7j for mixed SU(N) vertex model for the case N = 4. 

An exact solution can be found also in the non critical region which corresponds to 
fugacity n > 2. Since the fugacity is defined by n = 2 cos7r6', this region corresponds 
to imaginary values for the parameter 6. For convenience we assume 6 real so that the 
fugacity is now given hj n = 2 cosh 716', and redefine the function (|35|) and (|43|) by: 

I sinh^(c + xi) , , ^ 

Now in the thermodynamic limit, the roots corresponding to the largest eigenvalue are 
symmetrically distributed in the finite interval [—1/6', 1/6'] . This time the integral equation 
can be solved by Fourier series. The Fourier components are: 

p^(m, 6') = p^(m, ^^) = - cosh(m vr 6*) sech(2 m vr 6*) p^(m, 6') = - sech(2 m vr 6*) (58) 

The root density functions: 

Q +00 Q +00 

Pw{x,9) = - ^ p^(m, 6*) exp(im7r 6* x) pu,v{x,9) = - ^ p:;^(m, 6') exp(im7r 6' x) 

m=— 00 m=— 00 

(59) 

These can also be written in terms of Jacobi elliptic functions: 

p^(x, 6) = ^JIM. dn [x K{k) 6, k] (60) 

2 TT 

Pu,v{x,e) = — ^ {dn[x K{k) e + i K{k) e, k] + dn[x K{k) 9 ~ i K{k) 9, k]) (61) 

where the modulus k and the complete elliptic integral of the first kind K{k) are related 
to the phase 9 by: 

K{k) = r'\l - e sin^ 0]-V2 d(f) ^^i}'^'^ = 2 9 (62) 







K{k) 
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The free energy: 



. ^ c"'" ^ sinh(m VT ^) 1, ^1 l-g'" , , 

/oo W = V + E uo m = -7 log^ + E ^TT^ 63 

2 ^-^ mcosh{2m7T 6) 4 ^-^mq"^ + q ^ 

m=l ^ ' m=l 

where g = exp(— 2 tt ^^). Taylor expanding the summand around g = and summing with 
respect to m we get the following product expansion for the partition function: 



p=i 

For g ^ 1 we get the partition function for the four colorings model: 



U™oo^ ^ 11 (1 - g4p-3) (1 _ ^4p) 



11 2p(4p-3) v^,r3/2 ^^"^ 



in agreement with the expression (jSSI)- 

We have all the ingredients to discuss the singular behavior of the free energy. We can 
read the non analyticity of the partition function from expression (j64p . It appears that 
has a sohd wall of singularities on the unit circle ^] . Following ^Jj EH] we calculate 
the analytic continuation of the critical free energy into the non critical region so 
that we can extract the singularity of the free energy near the critical point n = 2. We 
find that near n = 2 the free energy has an essential singularity: 

fsing oc exp (66) 

This means that at the critical point n = 2 the model undergoes an infinite order phase 
transition of the Kosterlitz-Thouless type. The same type of transition has been found in 
[2ni for a loop model on the hexagonal lattice. Even if the free energy exhibits an essential 
singularity at n = 2, the right and left derivatives, of generic order, approach the same 
values, in the limit n — 2^. 

For the first derivative we get: 



d/ooH 



dn 



It has been shown that the first order derivative of the free energy with respect to the 
fugacity is related to the average loop length L [2]. We found L = 12, which is three times 
the minimal loop length allowed by the square lattice, in agreement with the conjecture 
made in 0. 
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6 Conclusions 



We have studied the bulk properties of the FPL^ model in the case in which the two loop 
fugacities are the same, both in the critical (n < 2) and non critical {n > 2) region. The 
model undergoes an infinite order phase transition, of the Kosterlitz-Thouless type, at 
the critical point n = 2. The Ansatz we make on the wave function works only when 
the vertex weights are the same uj = rj. This is because the amplitude that describes the 
scattering of particles belonging to different (odd and even) families, does not factorize 
into two particles scattering amplitudes when uj ^ rj. The details are in the Appendix. 
It will be interesting to find the Yang-Baxter relation on this solvable line and hopefully 
extend it to the whole phase diagram. On the solvable line the central charge and the 
scaling dimensions can be calculated exactly using the non linear integral equation (NLIE) 
method. 



7 Appendix 

Replacing 0^ in Eq- (|2S1) by Eq. (j221) the consistency requirement becomes: 



Zi-N a ZjZk 

2\ m+N-1 



(2\ m-i-jv— 1 ^ 







n=m+l 

In order to eliminate A, plug in the value of A ipm and A ipm+i given by the rhs of Eq. (fTUj) : 



J I I ^ ■ •'^ n A 

,4 



m+N-1 



Zi-]\f CL ZjZj^ Zi-]\[ 00 CL ^ 

Zm+l:n 



Zj{Zj^UJ + Zj) ^a„^„ +t^^/?n+l^n+l) "^^^^^ — (69) 

Zn+l:N+m-l 

n—m 
m+N 

/ -1 I 4\ / -2 I , 2 n I \ Zfn+2:n 
Zi [Z,^ + ZkUJ ) ^ [U anWn+i^ Pn+1 Vn+l) 



3 \ k 

^ -, Zn+l:N+r, 
n=m+l 

m+N-l 

-2,2 \ I -2 1^1 I -2 I , 2 o I \ Zm+2:n 



/ „2\ m+iv-i 
+ (r?"^ + ZjZk) Tj^'^Zk H V (t^~^ + ^n+l) 



We look for a sufficient condition on the scattering amplitude A j^... which ensures the 
validity of the equality (jU^ . We start by considering the terms, in the sums, that involve 
ipj for j 7^ m, m + 1 and require their single contribution to be zero. Working out we get 
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the following condition: 



1 1 2 

_ Zj 



(70) 



When 7] = uj the expression in parenthesis is antisymmetric and the sum vanishes if the 



scattering amplitude is symmetric A 



.jk. 



A 



.kj. 



Sorting out the terms ipm and ipm+i 



in (j69j) we got a more complicate expression which also goes to zero on the previous 
conditions, but not for the more general solution of (jTO)). 
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